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We study the combined impact of a colored environmental noise and demographic noise on the 
extinction risk of a long-lived and well-mixed isolated stochastic population which exhibits the 
Allee effect. The environmental noise modulates the population birth and death rates. Assuming 
that the Allee effect is strong, and the environmental noise is positively correlated and Gaussian, 
we derive a Fokker-Planck equation for the joint probability distribution of the population sizes 
and environmental fluctuations. In WKB approximation this equation reduces to an effective two- 
dimensional Hamiltonian mechanics, where the most likely path to extinction and the most likely 
environmental fluctuation are encoded in an instanton-like trajectory in the phase space. The 
mean time to extinction r is related to the mechanical action along this trajectory. We obtain 
new analytic results for short-correlated, long-correlated and relatively weak environmental noise. 
The population-size dependence of r changes from exponential for weak environmental noise to no 
dependence for strong noise, implying a greatly increased extinction risk. The theory is readily 
extendable to population switches between different metastable states, and to stochastic population 
explosion, due to a combined action of demographic and environmental noise. 
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I. INTRODUCTION 



A long-lived isolated stochastic population ultimately 
goes extinct via a large fluctuation: an unusual chain of 
deleterious events resulting from the demographic noise 
(the intrinsic discreteness of individuals and random na- 
ture of birth-death processes) and environmental varia- 
tions, see Ref. [1] for a recent review. It is important to 
understand how the interplay of environmental and de- 
mographic noises determines the mean time to extinction 
(MTE) 0. Early models postulated that the environ- 
mental noise, which modulates the birth and death rates 
of the population, is delta-correlated in time 0, 0] • Later 
on, population biologists realized, mostly via stochas- 
tic simulations, that temporal autocorrelation, or color, 
of environmental noise may have a considerable effect 
on population extinction [l], 0] ■ These insights inspired 
physicists who developed a theoretical framework for the 
analysis of a joint action of demographic and colored en- 
vironmental noise on extinction of an established popula- 
tion whose dynamics follows a simple stochastic logistic 
model 0. This theoretical framework provided a trans- 
parent way of evaluating the MTE and finding the op- 
timal environmental fluctuation that determines the op- 
timal (most likely) path of the population to extinction. 
The theory of 5] predicted the MTE in different regions 
of a two-dimensional "phase diagram" whose axes are the 
properly rescaled intensity (or, alternatively, variance), 
and the correlation time of the environmental noise. It 
tracked how the population-size dependence of the MTE 
changes from exponential with no environmental noise to 
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a power law for a short-correlated noise and to no depen- 
dence for long-correlated noise. It also established the 
validity domains of the white-noise limit and adiabatic 
limit. (In the adiabatic limit the environmental noise is 
assumed to vary very slowly compared with the relax- 
ation rate of the population toward the attracting fixed 
point of the deterministic rate equation.) 

The simple logistic models adopted in Refs. 0-01 do 
not account for the demographic Allee effect, by which 
population biologists mean a host of effects leading to 
an effective reduction in the per-capita growth rate at 
small population size ■ When the Allee effect is signifi- 
cant, a non-zero critical population size for establishment 
arises. If the initial population size is smaller than the 
critical size, the population quickly goes extinct. If the 
initial population size is greater than the critical one, 
an established population appears. Population biologists 
have argued that the Allee effect may influence, in a sig- 
nificant way, the population extinction risk due to the 
demographic and environmental noise Q- No satisfac- 
tory theoretical framework, however, has been developed 
until now. 

The present work attempts to close this gap. We for- 
mulate a minimal theoretical framework for this prob- 
lem by considering a simple set of stochastic reactions 
which mimics the Allee effect in a well-mixed popula- 
tion. The per-capita rates are modulated by a posi- 
tively correlated Gaussian noise with given magnitude 
and correlation time. We assume that the Allee effect is 
so strong, that the established population size, as pre- 
dicted by the deterministic rate equation, is close to the 
critical population size for establishment. In this limit 
(that is, close to the saddle-node bifurcation of the deter- 
ministic rate equation) a Fokker-Planck equation can be 
derived, which accurately describes the time evolution of 
the joint probability distribution of the population sizes 
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and environmental fluctuations. Throughout this work 
we assume that both the environmental noise and the 
demographic noise are weak, so that the MTE of the 
population is very long compared with the characteristic 
relaxation time predicted by the (noiseless) determinis- 
tic rate equation for this population. This enables us 
to use a small-noise approximation due to Freidlin and 
Wentzell essentially, a dissipative variant of WKB 
approximation. The WKB approximation reduces the 
Fokker-Planck equation to an effective two-dimensional 
classical mechanics. The optimal path of the population 
to extinction and the optimal environmental fluctuation 
are encoded in an instanton-like trajectory in the Hamil- 
tonian phase space of this classical mechanics, while the 
MTE is related to the mechanical action along the in- 
stanton. 

We solve the effective mechanical problem, and obtain 
analytic estimates for the MTE, perturbatively in three 
different limits: of short-correlated, long-correlated, and 
relatively weak environmental noise, for a population ex- 
hibiting a strong Allee effect. We also find, in each of 
these limits, the optimal (most likely) path of the popu- 
lation to extinction and the optimal environmental fluc- 
tuation. We complement our analytic results by solv- 
ing numerically the equations of motion of the effective 
classical mechanics. We find that the Allee effect has a 
strong impact on the MTE. It was discovered more than 
30 years ago by Leigh [3, |3] that, without the Allee ef- 
fect, a strong uncorrelated (white) environmental noise 
changes the population-size dependence of the MTE from 
an exponential to a power-law with a large exponent. We 
show here that, in the presence of the strong Allee effect, 
no power law appears. Here the population-size depen- 
dence of the MTE changes from exponential for weak 
environmental noise to no dependence for strong environ- 
mental noise. Our theory is readily extendable to popu- 
lation switches between different metastable states, and 
to noise- induced population explosion, due to a combined 
action of demographic and environmental noise. Where 
possible, we compare our results with previous ones. 

We reiterate that both demographic and environmen- 
tal noises are weak in our theory. Therefore, when we call 
the environmental noise weak or strong, we only mean 
that it is weak or strong compared with the demographic 
noise. 

Here is a plan of the remainder of the paper. Sections 
II and HI include preliminaries. In Section II we intro- 
duce a simple model of long-lived stochastic population 
which exhibits the Allee effect and ultimately goes extinct 
because of demographic noise. We start with the deter- 
ministic limit of the model and then outline its stochastic 
behavior, focusing on the limit of strong Allee effect. As 
a preliminary, we present in Section II a calculation of 
the MTE based on WKB approximation. In Section HI 
we add environmental noise to the model: first white 
noise, and then colored noise. In Section IV we evalu- 
ate the MTE of the population under the simultaneous 
action of environmental and demographic noises. Differ- 



ent subsections of Section IV deal with different limits: 
of short-correlated, long-correlated and (relatively) weak 
environmental noise. Section IV also includes a brief dis- 
cussion of (relatively) strong environmental noise where 
our results for the MTE coincide with previously known 
results. Our main findings are summarized in Section V. 

II. STOCHASTIC POPULATION WITH THE 
ALLEE EFFECT: PRELIMINARIES 

In the absence of environmental noise, a stochastic 
population exhibiting the Allee effect can be mimicked 
by three elementary reactions describing binary repro- 
duction, its inverse process and linear decay Q: 

2A^3A, A A 0, (1) 

with the (constant) reaction rates A, a and /j,. 

A. Deterministic Rate Equation 

The deterministic (or mean-field) theory only deals 
with the mean population size (the number of A's in the 
system), which is assumed to be large: n{t) ^ 1. The 
deterministic rate equation has the form 

f{n) = -fin+^n^ -^n^. (2) 

When — 1 — 8/iO'/(3A^) > 0, this equation has three 
fixed points and, therefore, describes a significant Allee 
effect. The fixed points uq ~ and — K{1 + 5) are 
attracting, the fixed point n_ = K{1 — S) is repelling. 
The parameter K = 3A/(2cr) plays the role of carrying 
capacity, as it sets the scale of the established population 
size. We will assume if 1 throughout this work. Equa- 
tion ^ describes an overdamped dynamics of a classi- 
cal particle with "coordinate" n in the effective potential 
U{n) = —J"'f{x)dx, see Fig. [TJ The fixed point 7i_ 
corresponds to the critical population size for establish- 
ment, whereas n+ corresponds to the established popu- 
lation. That is, according to the mean-field theory, once 
the initial population size exceeds n_, the population 
size will approach the fixed point tlj^ and remain there 
indefinitely. 

B. Stochastic Description and WKB 
Approximation 

The mean-field theory, however, disregards fluctua- 
tions of the population size around n = n+ . These fluctu- 
ations are caused by the demographic noise coming from 
the discreteness of particles and from the stochastic char- 
acter of the reactions H]). As K ^ 1, these fluctuations 
are typically small. However, a rare large fluctuation (an 
unusual chain of deleterious reactions) ultimately arises 
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FIG. 1: Effective potential U{n) = — /" f{x) dx correspond- 
ing to Eq. ((2| with 5^ > 0. The evolution of n{t) coincides 
with the motion of an overdamped particle in this potential. 
Above n = n_ the population gets established at n = n+, 
whereas below n — n_ the population goes extinct. 

and drives the population to extinction. Indeed, with 
the death of the last particle, there is no mechanism that 
would replenish the population. Because of the Allee 
effect, it is sufficient for the fluctuation to bring the pop- 
ulation below the critical point n = n_, whereupon the 
population goes extinct essentially deterministically 

Fluctuations of the population size are encoded in the 
probability Pn(t) to have, at time t, a population of n 
particles. Th e dy namics of Pn{t) is governed by the mas- 
ter equation p^lTl| 

Pn = HPn = Xn-lPn-1 + fJ-n+lPn+1 - + ^n)Pn, (3) 

where 

Xn(n — 1) , crnfn — — 2) 
A„ = — and ^„ = — ^ ^ + ^in (4) 

are the effective birth and death rates. After a short 
relaxation time, determined by Eq. a long-lived 

metastable distribution is formed where Pn{t) becomes 
sharply peaked at n = ri+ with an exponential decay to- 
wards n = n_ and an almost flat tail at < n < n_ 
^ . This almost flat tail determines the very slow "prob- 
ability leakage" into the absorbing state at n = 0. The 
leakage is described by the (exponentially small) lowest 
positive eigenvalue 1 /t of the operator H: 

P„(t) ~ 7r„e-*/^ n>0. (5) 

Here 7r„ is the lowest excited eigenstate of H: 

HlTn = (1/r) TTn (6) 

that can be identified with the quasi- stationary distribu- 
tion (QSD). In its turn, r is equal to the MTE, as 

PoW ^ 1 - e-*/^ (7) 

see Ref. for detail. 

Employing the large parameter K ^ 1, one can ac- 
curately calculate the MTE and QSD in this and many 



other one-population models @. Our present strategy, 
however, is different. We will make an additional assump- 
tion of a very strong Allee effect where the critical pop- 
ulation size n_ is relatively close to the established pop- 
ulation size n+ as described by the deterministic system 
Equivalently, the system is close to its saddle- 
node bifurcation corresponding to the appearance of the 
fixed points n± . In this limit a whole class of population 
models behaves in a universal way @. Furthermore, in 
this limit Pn(t) varies with n sufficiently slowly, and the 
van Kampen system size expansion (lOl [TT| becomes an 
accurate and controllable procedure. This procedure re- 
places the master equation ([3]) by a Fokker-Planck equa- 
tion, see below. Then the Langevin equation, equivalent 
to this Fokker-Planck equation, can be conveniently used 
for the introduction of environmental noise. 

In our example of three reactions, the strong- Allee- 
effect regime is achieved when 5 <C 1. In this case the 
MTE becomes ^ 

The approximate Fokker-Planck equation can be derived 
from Eq. ([3]) in a standard manner [l^, [lH • It reads 

d 1 d"^ 

Pn = r- [(A„ - /i„) Pn] + T,-r^ [{^n + M") Pn] , 

dn I dn^ 

where n ^ 1 is treated as a continuous variable. Rescal- 
ing time, t — crK'^t/Q, and the population size, q ~ n/K, 
we obtain 

V = {q [{q If 5'] V]' + ^{q [{q + if - 5^] V]" , 

(9) 

for the continuous probability distribution V{q,t). The 
overbar in t is omitted; the primes denote derivatives 
with respect to q. 

Since after a short transient V{q,t) becomes sharply 
peaked at q = 1 + 5, and (5 <C 1, we can simplify Eq. (|9]) 
by putting q = 1 everywhere except in the combination 
g — 1, and by neglecting in the second term on the 
right. The resulting equation is 

^ = {[ii-^r'S^]vy+^v". (10) 

This Fokker-Planck equation is equivalent, see e.g. 
Ref. [ll| , to the following Langevin equation: 

q = -{q-lf+S^ + ^r^4t), (11) 

where rjdit) is a white Gaussian noise with zero mean, 
{Vd{t)) = 0, and {rjd(ti)rjd(t2)) = S{ti - t2), whereas the 
subscript d stands for demographic. That is, close to the 
saddle-node bifurcation, the demographic noise is effec- 
tively Gaussian, white (that is, uncorrelated) and addi- 
tive. From now on, when discussing the correlation prop- 
erties of a noise, we will always mean the environmental 
noise. 
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Equations ^TU\\ and ([TT|) hold (up to rescaling, and 
close to the saddle-node bifurcation) for a whole class of 
single-population models which exhibit the AUee effect. 
Furthermore, these equations represent a truly paradig- 
matic model of escape from the vicinity of an attracting 
fixed point due to a weak additive white Gaussian noise. 
This model has appeared in numerous contexts, and the 
mean time to escape in this model is well known 10, 11]. 
We will proceed, however, as if we were unaware of these 
classical results. This is because we want, as a prelim- 
inary for the following material, to briefly outline how 
to evaluate the mean time to escape by using the weak- 
noise WKB approximation due to Freidlin and Wentzell 
see also Refs. [13, 113 ■ As we will see shortly, this 
approximation is readily extendable to the situation of 
our interest where weak demographic and environmental 
noises are both present. 

WKB approximation predicts the mean time to es- 
cape which coincides with the MTE. We look for the 
solution of Eq. (ITUl) as P{q,t) — •7r(g) exp(— t/r) and 
assume that t is exponentially large with respect to 
the parameter K, see Eq. ([5]). This justifies a quasi- 
stationary formulation for 7T{q). Then the WKB ansatz 
Tr{q) = exp[KS{q)] yields, in the leading order in 1/K, 
a stationary Hamilton- Jacobi equation H{q, dqS) ~ 0, 
where 

H{q,p)=2p'-[{q-ir-S^]p, (12) 

and p is the "momentum" canonically conjugate to the 
"coordinate" q. We should only deal with zero-energy 
trajectories. One type of zero-energy trajectories also 
have a zero momentum, p — 0. For p ^ the Hamil- 
ton equations read q = —{q — 1)'^ + 6'^ , p = 0, so these 
are (deterministic) relaxation trajectories. The escape is 
encoded by an activation trajectory: a zero-energy tra- 
jectory with p 7^ 0. This trajectory, 

p^Po{q) = {l/2)[{q-ir-5'] (13) 

is an instanton or, in a more mathematical language, a 
heteroclinic connection that exits the fixed point q = 
1 + S, p = and enters the fixed point q = 1 — S, p ~ 0. 
Then the population size flows toward q = along a 
deterministic trajectory which does not cost action, see 
Fig-in (The latter segment corresponds to the aforemen- 
tioned almost flat tail of the probability distribution.) 
Therefore, in the leading WKB order, we only need to 
calculate the mechanical action along the the activation 
trajectory from q = 1 + S to q = I — S: 

So= f \o{q)dq^ls\ (14) 
Jl+S 

As a result, the MTE due to the demographic noise is, 
up to a pre-exponent 

TO exp (I . (15) 



This result (which is of course well known) agrees with 
the more accurate result presented in Eq. ([5]). The 
quantity {2/3)K6^ is nothing but the Arrhenius factor 
AUq/Q, where Q = 2/K is the effective temperature [see 
Eqs. (Uni) and ^],AUo^Uo{q = l-S)-Uo{q^l + S) 
is the potential barrier height, and Uo{q) ~ {l/3){q — 
1)^ — S^q is the potential corresponding to the force 
/o(g) = -(g-l)2 + 52. 

Now we can see that the large parameter of the WKB 
theory is KS^ » 1. Actually, this could have been seen 
directly from Eq. (ITUl) : by rescaling (g — l)/6 q and 
6t t, one obtains the equation 

containing a single parameter KS^. When this parame- 
ter is large, the quasi-stationary distribution is sharply 
peaked around the attracting fixed point q = 1, thus val- 
idating the WKB approximation. 
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FIG. 2: (Color online) Zero-energy phase trajectories of the 
Hamiltonian (|12p . The daslied line shows the p = trajectory. 
The solid curve depicts po(<?) from Eq. (|13l) . The area of the 
shaded region is equal to So from Eq. (|14p . The thick line 
shows the path to extinction. 



III. INCORPORATING ENVIRONMENTAL 
NOISE 

Environmental variations affect the population dynam- 
ics by modulating the birth and death rates. As a re- 
sult, the bifurcation parameter 5^, which enters Eq. (jlip . 
becomes time-dependent: S^{t) ~ — (.(t), where the 
zero-mean random process ^(t) is independent of the de- 
mographic noise r]d{t). Now Eq. (fTTj) becomes 

9 = -(9-l)'+^o-e(0 + /|%W- (17) 

In fact, making the reaction rates noisy will in general 
affect not only S'^ but also K. However, if the environ- 
mental noise is sufficiently weak, an account of this effect 
only leads to a subleading correction that we will ignore. 

In contrast to the Langevin equation that describes a 
combined action of the demographic and environmental 
noise in the absence of the AUee effect [l[, the demo- 
graphic and environmental noises in Eq. (jl7p are both 
additive. 
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A. White Noise 

The nature of environmental noise manifests itself in 
the properties of the random process ^ {t) . The simplest 
environmental noise to consider is a white noise of a 
given intensity D: ^(t) = \/2Dr]e{t), with r]e{t) having 
the same properties as rjd{t)^ the two being independent. 
Now Eq. (HT]) reads 

q = -{q-lf + 51 + \I^Vd{t) + V2Drje{t). (18) 

As r]d and rje are statistically independent Gaussian pro- 
cesses, their weighted sum is another Gaussian, and we 
obtain 

q = -{q-lf + 6l+i^- + 2Dj rjtit), (19) 

where Tyt (i) is a white Gaussian noise with the same prop- 
erties as rid{t) and f?e(i)- Equation (|19p coincides with 
Eq. (fTT|) . except for a greater noise intensity. If the com- 
bined noise is still sufficiently weak, we can again use 
WKB approximation, or simply replace K in Eq. (jl5l) by 
K/{I+1), where 



/ -DK 
2 



(20) 



is the ratio of the environmental and demographic noise 
intensities. This corresponds to the MTE 



'^whitc 



with 5, 



white 



2^3 



3(/+l) 



(21) 



The reduction of the purely demographic action Sq by a 
factor of /-|- 1 has a great impact on the MTE, especially 
if / is large. When / >• 1, 5'whitc scales as l/I. As 
a result, the demographic parameter K cancels out in 
the expression for the MTE, and r ^ exp[(45^)/(3-D)] is 
dominated by the environmental noise. 



B. Colored Noise 

Let us now return to Eq. pT)) and choose £,{t) to be 
a colored (positively correlated) Gaussian noise, as mod- 
eled by the Ornstein-Uhlenbeck random process [ll| . The 
autocorrelation function of the Ornstein-Uhlenbeck ran- 
dom process is (^(^1)^(^2)) = (D/tc) exp{-\ti - i2|/Tc), 
where Tc is the correlation time of the noise, and D is the 
noise intensity. As Eq. (|17p is written in rescaled vari- 
ables, Tc and D are also assumed to be properly rescaled. 
The variance of the Ornstein-Uhlenbeck noise is equal 
to Dtc- The Ornstein-Uhlenbeck process can be conve- 
niently represented by a Langevin equation. 



When Tc tends to zero, Eq. (|22)) turns into ^{t) = 
\J2D 77e(i), and the white-noise limit is recovered. For fi- 
nite Tc we have to deal with two coupled scalar Langevin 
equations (I17p and (j22p with mutually independent white 
Gaussian noises ?/e(0 and rid(t). These Langevin equa- 
tions are equivalent, see e.g. [ll|, to a Fokker-Planck 
equation for the joint probability distribution V{q,^,t) 
of the rescaled population sizes q and the environmental 
fluctuations ^: 



dV d 

9f = 9^^[('^ 



^)'-^o]V}+^^ + -^ + 

-g-^m^^^- (23) 



IV. WKB ANALYSIS 

We assume throughout this paper that both envi- 
ronmental and demographic noises are sufficiently weak 
(we will obtain the corresponding criteria a posteriori). 
Mathematically, this means that the coefficients of the 
two diffusion terms in Eq. (1^^ are small. In this case 
the time history of V{q, ^, t), as described by Eq. (^51 . is 
the following. After a relatively short transient (typically 
of duration ^ l/(5o), V{q,^,t) develops a sharp peak at 
q = 1 + (5o, C = 0, as predicted by Eq. with if -> 00 
and D —i' 0. The small diffusion terms arrest the peak 
growth so that quasi-stationarity sets in. At the same 
time, the probability very slowly leaks toward the point 
q = 1 — Sq, ^ — beyond which the distribution is al- 
most flat, which corresponds to a quick escape toward 
the absorbing state at q = 0. This late-time dynamics is 
described by the eigenstate of the Fokker-Planck operator 
with the (exponentially small) lowest positive eigenvalue. 
Combining this knowledge with the Freidlin-Wentzel (or 
WKB) ansatz for the quasi-stationary distribution, we 
can write 

riq,U)^AQ,Oe'K ^(<Z,e)=e-^'^(''^«). (24) 

Then Eq. ([23]) yields, in the leading order in 1/K, a 
Hamilton- Jacobi equation with two degrees of freedom: 
H{q,£^,dqS,d(^S) ~ 0, with effective Hamiltonian 

H = 2p'- [{q If S^q] p-^p+^P^^ l^P. (25) 

Tc Tc 

The two momenta p and P are conjugate to the "coordi- 
nates" q and ^. The Hamilton equations are 

q = si^c~{q-ir + 4p, + 

Tc T,^ 



p^2{q-l)p, 



1 

P= —P + p. 



Tc 



■Ve{t). 



(22) 



Now we should look for an instanton: a zero-energy but 
non-zero momentum trajectory which exits, at t = —00, 
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the fixed point {q,^,p,P) = {1 + Sa,0, 0, 0) of tliis Hamil- 
tonian flow, and enters, at < = oo, tlie fixed point 
{q,^,p,P) ~ (1 — (5o,0,0, 0). Once tlie instanton is found, 
we can compute tlie action along it. 



S^S^ + S^^ J {pq + P^]dt, 
and evaluate the MTE from 



Let us apply one more rescaling transformation: 

q=iq-l)/So, e = e/^o, 
P^p/Sl P = P/So, 



t = dot. 



(26) 



(27) 



(28) 



Omitting the overbars in the equations of motion, we 
obtain 



9 = 






(29) 


P = 






(30) 


i = 




4/ 


(31) 


p = 


T 




(32) 



As we can see, the dynamics is controlled by two dimen- 
sionless parameters: / — DK/2 and T = SqTc, the latter 
being the ratio of the correlation time of the environmen- 
tal noise and the relaxation time of the system without 
noise. As we will see, sometimes it is more convenient to 
use the rescaled variance V = I/T of the environmental 
noise instead of the rescaled intensity /. 

Equations (|29p - ([5^ are Hamiltonian, as they stem 
from the Hamiltonian 

H = 2p^^{q^-l)p-^p+^P'-^^P (33) 

In the rescaled variables, the instanton connects the fixed 
points (1,0,0,0) and (—1,0,0,0). Denoting the action 
along this instanton by S, we can express the original 
action S that appears in Eq. as 5 = 6qS. As a 
result, Eq. ([27l) becomes 



(34) 



The two-dimensional Hamiltonian system p3p is in 
general non-integrable, as the only available integral of 
motion - the Hamiltonian itself - is insufficient for in- 
tegr ability. As a result, it is impossible to find the in- 
stanton analytically for arbitrary / and T. Perturbative 
solutions of different types are possible, however, in sev- 
eral regions of the (J, T) plane; seeking such solutions 
will be our main strategy for the remainder of the pa- 
per. We will refer to the limits of small and large T 
(with the criteria derived a posteriori) as the short- and 
long-correlated (environmental) noise, respectively. The 
limits of small and large / will be called the weak and 
strong (environmental) noise, respectively. 



Short-Correlated Noise 



Leading order in T <C 1 



For very small T the right-hand-sides of Eqs. ([31]) and 
include large factors. As a result, ^ and P quickly 
adjust to the current value of p{t) which slowly evolves 
with time: ^{t) ~ -AIp{t) and P{t) ~ -Tp{t). Plugging 
this ^{t) in Eq. we obtain 

q=l-q^ +A{I+l)p. 

This equation and Eq. (j30p are Hamilton equations, with 
the effective Hamiltonian 

H,iq,p)=2{I + l)p^ + {l-q^)p 

which coincides, up to rescaling, with Eq. p^ . The es- 
cape instanton satisfies 



Po{q) = 



g2 - 1 



2(/+l)' 

so the rescaled action in the leading order is 



So= 1^ ^"('^)^'^=3(/+l) 



(35) 



which, along with Eq. yields Eq. (gT]) for the MTE. 

One can also obtain explicit solutions for q{t) and p{t) 
[l5| in the leading order in T <C 1: 



qo{t) 



tanht, 



Poit) 



1 



2(J + 1) cosh^ t 



Correspondingly, 



m = 



21 



(/ + l)cosh^i 



Po{t) 



T 



2(/-f l)cosh^t 



(36) 



(37) 



where the subscript stands for the leading-order quanti- 
ties. As one can see, qo{t) does not depend on /, whereas 
the magnitude of the environmental fluctuation (t) goes 
up with / and then saturates: S,o{t,I ^ 1) = 2cosh~^t. 
The magnitudes of the momenta p and P go down as 
/ increases. This is expected on physical grounds: the 
stronger is the environmental noise, the smaller are the 
momenta needed for escape, leading to a smaller action 
and a shorter escape time. As both Po{t) and £,o{t) are 
even functions of time, the environmental noise does not 
contribute to the action in the first order of T. 

Figure [3] depicts qo,po,£,o and Pq versus time, 
with 6'^{t). Note that, for / > 1, the effective 
dependent bifurcation parameter S^{t) — Sl{l ~ 
becomes negative on a time interval around t = 
change of sign, however, occurs on the same time scale 
as that of the q- and p-dynamics, and does not lead to 
any qualitative change in the character of solution [l6j |. 



along 
time- 

This 
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FIG. 3: (Color online) Optimal trajectory (a) and opti- 
mal environmental fluctuation (b) for short-correlated noise. 
Po{t), £,o{t) and Po{t) are normalized by their extrema and 
denoted by overbars. Solid curves: qo{t) (a) and S,o{t) = 
(l/2)(l-hl/7)Co(t)andPoW = {2/T){I+l) Po{t) (b). Dashed 
curve: po{t) = 2(7-1- l)po(i) (a). Dot-dashed curve: 5^(t) (b). 



S. Subleading order in T <C 1 

Now we calculate the next-order correction to the 
white-noise result. Using T as a small parameter, we 
look for the solutions of Eqs. (|3T|) and (l32t as £,{t) — 
Ut) + T^Ut) + • • ■ and P{t) = TPi{t) + T^P2{t) + ... 
(the odd powers of T in the expansion of ^ turn out to be 
absent). We substitute these expressions into Eqs. ([51]) 
and (j32p and demand cancelation in every order in T. 
This procedure yields ^ and P expressed via p{t) : 



P{t) - -T[p(t)+Tp{t)+T^m + ---] 



(38) 
(39) 



where the leading-order terms coincide with those we ob- 
tained previously. Combining Eq. (p8)) with Eq. (p9)) . we 
obtain in the leading and subleading orders 



5 = l_92^4(/+l)(p + £p), 



(40) 



where e = IT"^ / {I + I) I. Equations ^ and (gOl) 
make a closed set and can be solved perturbatively in e, 



by setting q{t) = qo (t) + eqi {t) , pit) = po (t) + epi (t) with 
go and po from Eq. p6p . Eliminating pi we obtain, in 
the first order in e: 



91 



/ 6 



V cosh^ t 



4 Ui 



24 sinhi 
cosh^t ' 



We are looking for the forced solution of this linear equa- 
tion which obeys zero boundary conditions at t — ?> ±oo. 
This solution turns out to be elementary: 



4 sinh t 
cosh^ t ' 



(41) 



The corresponding forced solution for pi, that vanishes 
at t ^ ±cx), is 



2 sinh^ t 
(I + I) cosh'^ t' 



(42) 



Now we can calculate 5*^: the contribution of the {q,p) 
subsystem to the action (P5|: 



= I [Po<io + siPoQi +Piqo)] dt 



— oo 

2 



8/T2 



3(7+1) 15(7+1) 



Once p{t) is found up to the second order in £, the sub- 
leading corrections for ^ and P can be calculated from 
Eqs. ([55)) and ([5^ . respectively. We skip these formulas 
here and focus on calculating the important correction 
to the action coming from the (^, P) subsystem. Using 
Eqs. dMl) and we obtain 



= P£,dt^AIT 

J — OO 

i67r2 



po'dt + 0(T3) 



15(7+1) 



e'(T^). 



The total action is then 

S = Sq -\- — S'whitc 



47T^ 
5(7+1) 



(43) 



(44) 



where S^whitc = (2/3) (7 + 1)^^. That is, for an almost 
white noise, T < 1, the MTE ([34]) is longer, and so 
the extinction risk is lower, than for the white noise of 
the same intensity I. However, if we keep the variance 
constant, V = I /T = const, then T > reduces the 
MTE and increases the extinction risk, as follows from 
the leading-order result ([SSj): So = (2/3)(l + VT)'^ ~ 
(2/3)(l - VT). 

Figure |4| shows a comparison of the action from 
Eq. ([44l) with the results of our numerical calculations 
for 7 = 1 and different T. The numerical results were 
obtained by computing the instanton solution of the full 
set of equations (|29l) - ([5^ by a shooting method, and then 
evaluating the action integral in Eq. (j26p numerically, see 
Ref. 17,1 for details. 
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FIG. 4: (Color online) Analytic (lines) and numerical (circles) 
results for S/S^hnc versus T on a semi-logarithmic scale for 
I — 1. Left line: short-correlated noise theory, Eq. (|44|l . 
Right line: long-correlated noise theory, Eq. (|56|) . The inset 
shows, on a log-log scale, the small-T correction S/S^hitc — 1, 
see Eq. (|44p . The asymptote's slope is 2, so the behavior 
of the correction is evident. 



entering the forcing term, is given by Eq. (|47)) . but 
v = vit) = [1 — £,{1)]^/"^ is now time-dependent. How- 
ever, the time scale T, determined by the left-hand side 
of Eq. (|48l) , is supposedly much longer than the time scale 
of the forcing. Therefore, the forcing pulse, see Eq. (|47)) 
for p, can be approximated by a delta-function with the 
proper amplitude: 



m 



y2 y3 

Here we have replaced iy{t) by 



Sit). 



^0 



it = 0) = ^i-at = o); 



(49) 



(50) 



the corresponding criterion will appear shortly. The so- 
kition of Eq. (gS]) is 



( 

T 



\t\/T 



(51) 



In its turn. 



B. Long-Correlated Noise 



P{t) = 



v^e^l^ t < 0, 
t>0. 



(52) 



Here analytic progress is possible due to time-scale sep- 
aration. Indeed, for sufficiently large T (we will obtain 
the criterion a posteriori) the right-hand-side of Eq. (j3ip 
is small. Therefore, £,{t) varies slowly (adiabatically) 
compared with q{t) and p{t), and can be treated as con- 
stant when dealing with the fast (q,p) sub-system. Equa- 
tions ([^ and ([5(1)) with ^ — const are Hamilton equa- 
tions with the Hamiltonian 



Ho ^ 2p^ - {q^ - ly^) p, 



(45) 



where we have defined = 1 — ^. This Hamiltonian co- 
incides with that of Eq. ([T^ up to rescaling. We imme- 
diately obtain the instanton solution in parametric form: 



The time-dependent solutions are 

q{t) = —i/tanhiyt, 

1 



Pit) 



(46) 



(47) 



2 cosh lyt 



The characteristic fast time scale is l/i', with a yet un- 
known ly. We assume here (and will check a posteriori) 
that = 1 - ^ > 0. 

Now we turn to the slow sub-system {£,,P)- Differen- 
tiating Eq. (|31l) with respect to time and using Eq. (|32]) . 
we obtain an exact linear second-order equation for (t) : 



m-^-^-^Pit), 



(48) 



which has to be solved with the boundary conditions 
^(i — > ±oo) = 0. In our adiabatic approximation p{t), 



What is left is to find vq by equating £^{t 
Eqs. dSni) and dHl)- We obtain 



1 - V. 



0) from 



(53) 



As one can check, ma.x£^{t) ~ ^(0) < 1 with this fo, as we 

assumed. Figure [5] shows the analytic results for q,p,£ 
and P versus time, along with S'^(t). The same figure 
also shows our numerical results for the same / and T. 
The corner singularity in ^(t) and the jump in P{t), both 
observed at t = 0, are the approximation price we have to 
pay for replacing p{t) by the delta-function in the forcing 
term of Eq. (|48|. These singularities do not cause any 
problem in the action calculations that we now present. 

With Eqs. dSU and ([52]), the calculation of is 
straightforward: 



Pidt=^ = V [^VV^ + 1-V 



(54) 



The calculation of Sq, with q and p from Eq. (1471) . sim- 
plifies once we notice that the integral of pq over time 
is mostly gathered in a narrow time interval of width 
~ l/i/Q around t = 0. Within this interval one can re- 
place i'{t) by vq - the same replacement as in Eq. (H^ - 
and obtain 



5,= 



pqdt^^iy^ = '^(Vv^ 



1 - V 



The total action is 
S = S, 



(55) 



(56) 
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time l/t'o is short compared with the slow time T. Using 
Eq. (j53p . one can reduce this strong inequahty to 

T> max(l,\/7). (57) 

This condition, however, is insufficient. One also needs 
to demand that the variation of ^(i) during the fast 
time l/t'o be small compared with each of the terms of 
Eq. ((29|) : for example, with . This criterion can be 
written as <^(0)/i/o < ■ In view of Eqs. (gT]) and ((5T|) 
this criterion demands / ^ {vqT)^ which, after some al- 
gebra, boils down to T 3> max (J^/^, J^^^). [As one can 
check, the same criterion is required for the replacements 
of by vq in Eqs. and (|55p .] Combining this con- 
dition with Eq. (1571) , we obtain the adiabaticity criterion 
for the environmental noise: 

r>niax(l,/3/4) (58) 

or, in terms of the rescaled variance, 

r>max(l,y3)^ (59) 

For sufficiently large T Eq. (|56p agrees well with our nu- 
merical results, see the right solid line on Fig. |4l 

C. Weak Noise 

For sufficiently small / the problem can be solved per- 
turbatively. Let us split the Hamiltonian (j33p into un- 
perturbed and perturbed parts: H ~ Hq + I Hi, where 

Ho = 2p2 -{q^-l)p-^p- i^P, (60) 

Hi = (61) 

and / serves as the small parameter. Correspondingly, 
S = Sq + AS, where the small correction AS is propor- 
tional to /. 



FIG. 5: (Color online) Optimal trajectory (a) and optimal 
environmental fluctuation (b,c) for long-correlated noise for 
T = 16 and 1 = 4 (so that vq — 0.78). (a) Analytic results for 
q/uo (solid line) and p/i^o (dashed line) versus the fast time uot 
(main panel) and slow time t/T (inset). Numerical results are 
indistinguishable from analytic (not shown), (b) Numerical 
(solid line) and analytic (dashed line) results for ^/vo versus 
the slow (main panel) and fast (inset) times. Also shown is 
the analytic bifurcation parameter 5^ (t/T). (c) Numerical 
(solid line) and analytic (dashed line) results for P{t/T)/vo- 

Equation ([M]) with this S yields the MTE in this limit, 
up to a pre-exponential factor. 

What is the validity domain of the adiabatic approxi- 
mation that we have used? An obvious condition is the 
strong inequality vqT <^ 1 which guarantees that the fast 



1. Zeroth order 

The unperturbed, or zeroth-order problem is described 
by the Hamiltonian Hq, that is by Eqs. d^Hl-lISl]) with 
7 = 0. The zeroth-order equation for ^ is ^ = —S^/T; 
its only acceptable solution is ^ = = 0: no environ- 
mental noise. As a result, the zeroth-order equations ((29)) 
and (I30|) for q and p coincide with those without environ- 
mental noise, and their solutions, obeying the boundary 
conditions at t = ±oo, are 

qo{t) = - tanhi, po{t) = — 

2 cosh t 

The action, contributed by the {q,p) subsystem is Sq = 
2/3. Interestingly, the momentum Po(Oi conjugate to 



10 



(,Q{t) = 0, has a non-trivial behavior. It is described by 
the equation 

P-ip = L_ 

T 2cosh2t 
whose solution, vanishing at t — ±00, is 

t/T poo x/T 

Po{t) = / -—^dx. (62) 
^ Jt cosh X 

As ^o{t) — 0, this "ghost solution" does not contribute 
to the action, and = 2/3, coming from the (qo,Po) 
subsystem, yields Eq. We will need the "ghost so- 

lution" , however, in the first order calculations that we 
now present. 

2. First order 




0.5 5 50 



X 

FIG. 6: (Color online) The function <I>(a;) and its small- 
X (dashed) and large- a; (dot-dashed) asymptotics, 2/3 — 
(8/15)2;^ and 1/x, respectively, shown on a semi-logarithmic 
scale. 



The first-order correction to the action, AS', can be 
found by integrating IHi over the unperturbed trajecto- 
ries [5, 18, 19] 

/OO 
Hi[qo{t),po{t),^„{t),Po{t)]dt. 
-00 

Using Eqs. (|6T|) and (|62|) . we arrive at 
AS* 



I 

'2T2 



dt 



dx- 



-x/T 



dy- 



-V/T 



cosh X Jt cosh y 



-• (63) 



Evaluating this triple integral (see the Appendix for de- 
tails), we obtain 



AS- 



-I^T), 



(64) 



where 



$(x) = 



2x 



-x-l 



^p{x) — d^ InT [x) / dx^ , and r(a;) is the gamma-function. 
</? is the so called trigamma function: a special case of the 
polygamma function [l^l- The function <^{x) is plotted 
on Fig. ini along with its small- and large- a; asymptotics. 
Altogether, our weak-noise result for the action is 



3 r2 



— (—\-T 



(65) 



Using the small- and large- a; asymptotics of 1^3(2;), we 
can obtain simple formulas for S for short- and long- 
correlated noise 




T< 1, 
T> 1. 



(66) 



As one can easily check, the T ^ 1 asymptotic in Eq. ((66)) 
coincides with the / ^ 1 asymptotic of Eq. (j44|) obtained 
for the short-correlated noise. In its turn, the T ^ 1 
asymptotic in Eq. (I66|) coincides with the <C 1 asymp- 
totic of Eq. ([56l) obtained for the long-correlated noise. 
Figure [7] shows a comparison of Eq. (j65p with our numer- 
ical results for T = 1. One can see good agreement for 
sufficiently small /. 
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FIG. 7: (Color online) Analytic (lines) and numerical (circles) 
results for the action S versus / for T = 1. Left line: the 
weak-noise theory, Eq. (|65|) [also shown in inset (a)]. For 
/ ^ max(l, T) S must behave as S — g(T)/I, see subsection 
IIVDI The function g{T) is only known analytically for T <^ 1 
and for 1 < I^^^ < T < J. The right line is S = 0.9/1, that 
is g{T = 1) ~ 0.9. Inset (b): S versus I for very large I is 
displayed on a log-log scale. The slope of the asymptotic is 
— 1 as expected. 



Now we can determine the validity domain of the weak- 
noise approximation by demanding the strong inequality 
AS < So, or simply AS < 1. For T < 1 we have AS" ~ 
/, whereas for T > 1 we obtain AS ^ I/T. Therefore, 
the weak- noise approximation holds when / <^ max (T, 1) 
or, in terms of the rescaled variance, V <C max (l/T, 1). 
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D. Strong noise 

For strong environmental noise we only have partial re- 
sults which, as we will see shortly, are not new. How does 
S depend on / for very large /? Here the demographic 
noise becomes negligible compared with the environmen- 
tal noise. Therefore, the parameter K must drop from 
the exponent KS of the MTE in_Eq. This can 

only happen if S, and therefore S = S/6q behaves as 
S = g{T)/I. We can extract g{T) from our analytic re- 
sults in two different domains. For the short-correlated 
noise, T ^ 1, we can expand Eq. (|44|) at / 1. For the 
long-correlated noise, we can use the large- 1^ asymptotic 
of Eq. ([56]) . These procedures yield 




/> 1, T< 1, 
1 < < T < /. 



(67) 



These asymptotics can be compared with those obtained 
in 1989 by Bray and McKane [21|. They investigated es- 
cape of an overdamped particle from a smooth potential 
well U{x) solely due to an (Ornstein-Uhlenbeck) extrinsic 
noise with correlation time Tc and intensity D. The ab- 
sence of intrinsic noise in their setting corresponds to the 
limit of strong environmental noise in ours. Bray and 
McKane [2l| presented their result for the mean time 
to escape as Inr ~ s/D. To go over to our notation, 
we use Eq. to express D — 21 /K. As a result, 
Inr ~ Ks/{21), and our g{T) is related to their s as 
g{T) = s/(2c53). 

For short-correlated noise Bray and McKane arrived 

at 

pb 

s ^ U{a)-U{b) + T^ / dx[U"{x)]'^U'{x) 

J a 

l-b 

dx [U"'{x)f [U'{x)f + 0(t6), (68) 

where the fixed points a and b correspond to our 1 + 5 
and 1 — 5, respectively. Putting U{q) = Uo{q) = {q — 
1)^ /3 — Sgq, limiting ourselves only to the leading correc- 
tion 0{t^), and evaluating the integral in Eq. (|68)). we 
obtain 



16T^ 



4 

s = I - H 

,3 15 



where T = SqTc- This yields the first line in our E!q. (|67p . 

For long-correlated noise Bray and McKane [2l| ob- 
tained 



2 



[u'{d)Y 



(69) 



where d is the inflection point of the potential U{q), lo- 
cated between the points a and b. For our Uo{q) one has 
d = 1, and Eq. §9^ yields s = which leads to the 

second line in our Eq. (p7|) . 

For T ~ 1 analytic progress is difficult, as was already 
noticed in Ref. [2lj. Still, g{T) can be found numerically. 



see also Ref. [2]]]. For example, we found that g{l) ~ 0.9, 
see Fig. [T) Finally, the strong-noise limit corresponds to 
/> max(l,T), or F > max(l/r, 1). 



E. Phase diagram 

Table U summarizes our main analytic results for S ~ 
{KSq)~^ luT in different regions of the parameter plane 
{I,T). The regions themselves make a "phase diagram" 
which is shown in Fig. |S] in the (/, T) and (V, T) repre- 
sentations. 
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Adiabatic 
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1 - TfT i^fiw') - T - 1 





TABLE I: Action S in different parameter regions 
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FIG. 8: (Color online) Phase diagram of the system in the 
(/, r) (a) and {V, T) (b) representations. 

The validity domains of our results for the purpose of 
evaluation of the MTE can be found in each particular 
case by demanding that KS^S ^1. In the cases where 
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we calculated a small sub-leading term Si, a more strin- 
gent condition KSqSi ^ 1 is required. These criteria 
always hold for sufficiently large K and small D. 



V. SUMMARY 

We have evaluated the mean time to extinction (MTE) 
of a long-lived and well-mixed isolated population caused 
by an interplay of colored environmental noise and (ef- 
fectively white) demographic noise. We assumed that 
the population exhibits a strong AUee effect. We 
have obtained analytic results in the limits of short- 
correlated, long-correlated and (relatively) weak environ- 
mental noise, see Table 1. We have also established the 
validity domains of white noise and adiabatic noise on 
the parameter plane {I,T). As in the absence of the 
Allee effect, even a relatively weak environmental noise 
leads to an exponentially large reduction in the MTE. 
For a relatively strong environmental noise, this effect 
becomes dramatic. We have found, in the different lim- 
its, the most likely path of the population to extinction 
and the optimal environmental fluctuation (OEF) that 
mostly contributes to this path. 

For a relatively strong and short-correlated environ- 
mental noise the OEF temporarily changes the sign of 
the difference between the birth and death rates of the 
population. For long-correlated noise the OEF is such 
that this difference remains positive at all times, for any 
noise intensity. 

This theory is immediately extendable, close to the 
saddle-node bifurcation, to population transitions, due 
to a combined action of demographic and environmental 
noise, in two additional settings. The first setting is pop- 
ulation explosion. The second is population switches be- 
tween two non-empty states each of which, at the deter- 
ministic level, is linearly stable. Without environmental 
noise, these problems were considered in Refs. [3, [22| - 

Finally, it would be overly optimistic to hope that an 
analysis of a simple model that we presented here will 
resolve the long-time debate in population biology on 
"whether and under which conditions red noise increases 
or decreases extinction risk compared with uncorrelated 
(white) noise" Still, we believe that this analysis is 
a step toward resolving this debate. 
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APPENDIX: CALCULATION OF THE 
WEAK-NOISE INTEGRAL 

Here we present some details of the calculation of the 
triple integral in Eq. (1631) . Let us change the integra- 
tion order by moving the integral over t to the innermost 
position. Since the integration is over {(t,x,y) : x > 
t,y > t}, the X- and y-integration domains are (— oo, oo), 
whereas for any x, y the i-integration is from — cxi to 
u = Tjmi{x,y). The i-integration yields (T /2) e^^^^ . 
Now we split the integration domain of y into two sub- 
domains: y < X (where u = y) and to y > x (where 
u = x). We obtain 



AS^-— / dx 2- 

4-i J-oo cosh x 



dy- 



/-oo cosh^y 

roo ^(2x-y)/T 



cosh^ y 

27" J-oc J -oo cosh^ X cosh^ y 

' — X. Omitting the overbars and 
in order, we obtain 



Next we shift y: y = y — x. Omitting the 
changing the integration order, we obtain 





r 


2T j 


— oo 








— OO 


I 


r 


2r3 




/ 


1 


2^2 


2T 



pV/T 



fOC 



dx 

cosh^ X cosh^ {y + x) 



I —oo 

^ y cosh y 



, sinh'^ y 



1 ^ 
sinh^ y / 



— -\\- ' 

2T ^ j ^ (1 - 2T)2 



T- 1 



where 93(2) = lnr(z)/(iz^ is the trigamma function, 
and we have used the identity Lp[z — 1) — (z— 1)^^ ~ '('(z) 
a. 
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